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We conjecture that the confinement- deconfinement phase transition in QCD at 
large number of colors A'^ and Nf<^NatTj^O and ^ is triggered by the 
drastic change in 8 behavior. The conjecture is motivated by the holographic 
model of QCD where confinement -deconfinement phase transition indeed hap- 
pens precisely at T = Tc where 9 dependence experiences a sudden change in 
behavior. The conjecture is also supported by quantum field theory arguments 
when the instanton calculations (which trigger the 8 dependence) are under 
complete theoretical control for T > Tc, suddenly break down immediately be- 
low T <Tc with sharp changes in the 6 dependence. Finally, the conjecture is 
supported by a number of numerical lattice results. We employ this conjecture 
to study confinement -deconfinement phase transition of hot and dense QCD 
in large N limit by analyzing the 8 dependence. We estimate the critical values 
for Tc and /ic where the phase transition happens by approaching the critical 
values from the hot and/or dense regions where the instanton calculations are 
under complete theoretical control. We also describe some defects of various 
codimensions within a holographic model of QCD by focusing on their role 
around the phase transition point. 

Keywords: phase transition, large N expansion, 6 dependence, instantons. 



1. Introduction 

This talk is based on series of recent papers. Understanding the phase 
diagram at nonzero external parameters T, fj, is one of the most difficult 
problem in QCD. Obviously, this area is a prerogative of numerical lattice 
computations. However, some insights about the basic features of the phase 
diagram may be inferred by using some analytical approaches. In particular, 
some qualitative questions can be formulated and answered by considering a 
theory with large number of colors N. We conjecture that the confinement- 
deconfinement phase transition in QCD at large number of colors N and 



August 25, 2009 1:49 



WSPC - Proceedings Trim Size: 9in x 6in MN-proceedings 



2 

Nf N utT and /i 7^ is triggered by the drastic change in 9 behavior. 
This criteria is motivated by the observation that in holographic model of 
QCD the confinement -deconfinement phase transition happens precisely at 
the value of temperature T = T,. where 6 dependence experiences a sudden 
change in behavior. ^'^ Secondly, the proposal is supported by the numerical 
lattice results, see e.g. review article,'' which unambiguously suggest that 
the topological fluctuations are strongly suppressed in deconfined phase, 
and this suppression becomes more severe with increasing A'^. Finally, our 
new criteria is based on a physical picture which can be shortly summarized 
as follows. 

For sufficiently high temperatures T > Tc the instanton gas is dilute 
with density ^ g-iC^)^ which implies a strong suppression of the topo- 
logical fluctuations at large N where 7(T) > 0, see below for details on 
structure of j{T)— function. The calculations in this region are under com- 
plete theoretical control and the vacuum energy has a nice analytic behav- 
ior ~ cos ^e"'''^"^''^ as function of 9. At the critical value of temperature, 
T = Tc where j{T) changes the sign, the instanton expansion breaks down 
and one should naturally expect that at T = Tc there should be a sharp 
transition in 9 behavior as simple formula ^ cos 9 can only be valid when 
the instanton gas is dilute and semiclassical calculations are justified which 
is obviously not the case for T <Tc. Therefore, it is naturally to associate 
sharp changes in 9 behavior with confinement-deconfinement transition, 
just as in the holographic model. ^ 

The plan of the paper is as follows. I start in Section 2 by considering the 
phase transition in hot QCD.^ In section 3 the same technique is applied to 
the dense QCD. I argue that the confinement- deconfinement phase transi- 
tion happens at very large quark chemical potential fie ~ VNAqcd, where 
M = ^J-b/N is already properly scaled quark chemical potential.^ Finally, in 
section 4 some defects of various codimensions are described within a holo- 
graphic model of QCD by focusing on their role around the phase transition 
point."' 



2. Deconfinement Phase Transition in hot QCD at large AT. 

We start with analysis of nonzero temperature"'^ when a sharp transition in 
9 dependence at the phase transition is indeed observed in the holographic 
model of QCD. From quantum field theory viewpoint such a transition can 
be understood as follows. Instanton calculations are imder complete theoret- 
ical control in the region T > Tc as the instanton density is parametrically 
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suppressed at large N in deconfined region/ 

^nstW^e^^^cosf?, 7^ [iilnf-^) -1.861. (1) 

It is assumed that a higher order corrections may change the numerical coef- 
ficients in 7(r), but they do not change the structure of eq. H]). The critical 
temperature is determined by condition 7 = where exponentially small 
expansion parameter e~'^^ suddenly blows up and becomes exponentially 
large. Numerically, it happens at 



7 ^ 



3 (aqcd 



^ T,(A^-oo) ~0.53Aqci?, (2) 



where Aqcd is defined in the Pauli -Villars scheme. Our computations 
are carried out in the regime where the instanton density ^ exp(— 7iV) 
is parametrically suppressed at any small but finite 7(T) = e > when 
N = 00. From eq. ([1]) one can obtain the following expression for instanton 
density in vicinity of T > Tc, 

yi„,tW ~ COS0 • e-"^(^), 1 » » l/N. (3) 

where a = ^ and Tc{N — 00) ~ 0.53Aqc_d are estimated at one loop level. 
Such a behavior does imply that the dilute gas approximation is justified 
even in close vicinity of Tc as long as ^^^^k ^ Therefore, the 9 depen- 
dence, which is sensitive to the topological fluctuations is determined by 
([3]) all the way down to the temperatures very close to the phase transition 
point from above, T = Tc + 0{1/N). The topological susceptibility is order 
of one for T < Tc in confined phase while it vanishes ~ e~^^ — > for 
T > Tc in deconfined phase. Non topological quantum fluctuations on the 
other hand could be quite large in this region, but they do not effect the 
structure of eq. ([3]). 

There are three basic reasons for a generic structure ()f I2I3|1 to emerge: 

a. The presence of the exponentially large "T— independent" contribution 
( e.g. e+^-*^^ in eq. ([T])). This term basically describes the entropy of the 
configuration. It is due to a number of contributions such as a number of 
embedding SU{2) into SU{N) etc; 

b. The presence of the "T— dependent" contribution to V^nst(^') which 
comes from J n{p)dp integration. It is proportional to the standard factor. 



Aqcd\ '''' _ r 11.. . / ttT 
— — — — 



ttT J 



exp 



-— A^-ln 



3 V A 



^QCD 



(4) 
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c. The fermion related contributions such as a chiral condensate, diquark 
condensate or non-vanishing mass term enter the instanton density as fol- 
lows - {^^)^! - gAf-(«:in|{v;v>l)^ For K = ^ ^ this term obviously leads 
to a sub leading effects \/N in comparison with two main terms in the 
exponent ([T|) and will be neglected for the rest of this paper. 

The crucial element in this analysis is that both leading contributions 
(items 1 and 2 above) have exponential dependence, and therefore at 
iV — > oo for T > Tc the instanton gas is dilute with density e~^^ , 7 > 
which ensures a nice cos 6 dependence ([3]), while for T <Tc the expansion 
breaks down, and 6 dependence must sharply change at T < Tc. We have 
identified such sharp changes with first order phase transition. 

One more comment on this proposal. Our conjecture (that the 
confinement- deconfinement phase transition in QCD is triggered by the 
drastic change in 9 at the same point T = Tc) implicitly implies that the 
configurations which are responsible for sharp 9 changes must also dras- 
tically change their properties at T = Tc- For T > Tc the objects which 
describe the 9 behavior are the instantons, while at T < Tc we believe, they 
are the instanton quarks, the quantum objects with fractional topological 
charges zLl/N. Therefore, we interpret this transition as dissociation of the 
instanton into instanton-quarks which become the dominant quasi-particles 
at T < Tc, see more comments on this interpretation in concluding section. 

3. Phase Transition in dense QCD at large N. 

In this section we estimate the value of /Xc where the instanton expansion 
breaks down and therefore, the 9 dependence should experience a sharp 
change. According to our conjecture we should identify this position with 
the phase transition point. Similar arguments have been put forward pre- 
viously^ for numerical estimation of /ic for small N, Nf — 2, 3. Our goal 
here is quite different: we want to understand an analytical dependence 
of iJ,c{N,Nf) as a function of N,Nf at very large N and finite Nf ^ N 
in order to compare with results of ref.'' where the authors presented a 
very strong argument suggesting a very large fic ^ VN where the phase 
transition could happen. 

We follow the same logic as before, and study the 9 dependence in order 
to make a prediction about the phase transition point fic- In the regime 
H > fic the 9 dependence is determined by the dilute instanton gas approxi- 
mation. We expect that the expansion breaks down only in close vicinity of 
jjc at large N as it happens in our previous analysis with phase transition 
at T = Tc- According to the conjecture this point will be identified with 
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confinement- deconfinement phase transition point /ic- In the present case 
of analyzing fj,c ratfier tfian discussed previously we do not have any 
support from the lattice computations, nor from holographic models. Still, 
the basic governing principle remains the same. Therefore we identify the 
point where instanton expansion breaks down (and correspondingly a point 
where a simple cos6' sharply changes to something else) with the point /Xc 
where the phase transition happens. In our estimates below we assume that 
the color superconducting phase is realized in deconfined phase for all N, 
see e.g. recent review.® 

As we shall see below, the instanton density in deconfined phase has the 
following generic behavior, ~ cos 6* cxp [— A^7(/j,)], where 7(/i) ~ const. + 
0{1/N) in large N limit. Such a behavior implies that for any small (but 
finite) positive 7 > the instanton density is exponentially suppressed and 
our calculations are under complete theoretical control. In contrast: at ar- 
bitrary small and negative 7 < the instanton expansion obviously breaks 
down, theoretical control is lost as an exponential growth ~ exp (|7|iV) 
for the instanton density makes no sense. The 6 behavior must drastically 
change at this point. Therefore, the value of is determined by the fol- 
lowing condition, 

7(/U = yUc) = = cAqcd- (5) 

Our goal is to compute the coefficient c by approaching the critical point 
/ic from deconfined side of phase boundary. Therefore, we will be interested 
in the instanton density in the dilute gas regime at /i > /ic where analytical 
instanton calculations are under control. 

As is well known, the 9 dependence goes away in full QCD in both 
phases: confined as well as deconfined in the presence of the massless chiral 
fermions. However we are interested in the magnitude of the instanton 
contribution ~ T^nst {9) in deconfined phase rather than in dependence of 
full QCD. Precisely this coefficient triggers the point where the instanton 
expansion suddenly blows up. The sharp changes in T^nst(^) we identify 
with complete reconstruction of the ground state, drastic changes of the 
relevant gluon configurations, and finally, with confinement- deconfinement 
phase transition. To avoid identical vanishing of Finst (^) in the presence of 
massless fermions one can assume a non zero chiral condensate in deconfined 
phase or one can assume a non-vanishing masses mq 7^ for the fermions, 
or non- vanishing diquark condensate (ipfp) ^ to avoid identical vanishing 
of Vinst(^) for dense matter at large /i. None of these assumptions effects 
any numerical estimates given below in the limit N ^ 00, Nf <^ N, as 
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explained in the previous section. 

The result of the computations can be represented as follows, 

Must (0) ^e-^^ cos 0, ^=\^\n{^f-]-l.l\, ^? = Nii\ (6) 

where we introduced reduced chemical potential /2 = h/^/N and neglected 
all powers N''' in front of e~'^^ . The crucial difference in comparison with 
similar computation at nonzero temperature ([T]) is emerging of parameter 
jx instead of the original quark chemical potential = y/Nfi. It implies 
that the critical chemical potential where 7 changes the sign (and therefore 
where the phase transition is expected) is parametrically large fj,c ~ 
because /I^ ~ 1, see below for numerical estimates. The origin for this 
phenomenon can be understood from the following observation: the tem- 
perature -dependent factor in the instanton density before integrating over 
p is proportional to ^ iV while chemical potential enters this expression 
with factor ~ Nf ^ N, see e.g. the review paper^ . Therefore, a very 
large chemical potential /i ~ ^/NAqcd is required in order to achieve the 
same effect as temperature T Aqcd- The physics of this phenomenon 
can be explained as follows: at T ^ 1 a large number of gluons ^ N"^ 
can get excited while at ~ 1 only a relatively small number of quarks 
in fundamental representation ~ N can get excited. Therefore, it requires 
a very large chemical potential ji^ ^ N ui order for fundamental quarks 
play the same role as gluons do at T ^ 1. As explained above, the critical 
chemical potential is determined by condition 7 = where exponentially 
small expansion parameter e~''^ at /i > /Xc suddenly blows up at /i < /Xc- 
Numerically, it happens at 




1.1 



^ /ic(A^ 00) ~ 1.4 • Aqcd 




where ^.qcd is defined in the Pauli -Villars scheme. The topological sus- 
ceptibility vanishes ~ e~^^ — *■ for /i > /ic while it must be drastically 
different for /x < /x^ as 6* dependence must experience some drastic changes 
in this region as the instanton expansion breaks down, and therefore simple 
cos 9 dependence must be replaced by something else. It is very likely that 
the standard Witten's arguments (valid for the confined phase) still hold 
in this region ^ < in which case the topological susceptibility is order of 
one. 

The Aqcd in the Pauli -Villars scheme which enters our formula ([7]) 
is not well-known numerically. Therefore, for numerical estimates one can 
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trade Aqcd in favor oiTc{N = oo) at /i = estimated above (H]). Therefore, 
our final numerical estimate for fic{N = oo) can be presented as follows, 



Aic(iV = oo) ~ 2.6 • . /— -TcCA^ = (X),Ai = 0), iV/ < iV. (8) 



If one uses the numerical value for Tc{N = 3) ~ 260 MeV,^°'^^ one arrives 
to ndN = oo) ~ 690^y N/Nf MeV which is our final numerical estimate 
for the critical chemical potential where deconfincd phase transition is pre- 
dicted for very large N. Few remarks are in order: 

a. The most important result of the present studies is the observation that 
the confinement- deconfinement phase transition according to ([7]) happens 
at very large fic ^ Vn if Nf <C N. This is consistent with the results of'' 
where parametrically large scale for /ic ^ VN had been predicted. However, 
the technique of ref.^ does not allow to answer the question whether the 
transition would be the first order or it would be a crossover. Within our 
framework at ^ 1 and Nf <C N the entire phase transition line (which 
starts at T = Tc ~ Aqcd at /i = and ends at = /Xc ^ \/NAqcd at 
T = 0) is predicted to be the first order phase transition at large N and 
Nf <C N. This is because the nature for the phase transition along the 
entire line is one and the same: it is drastic changes of 9 dependence when 
the phase transition line is crossed. 

b. Our computations are carried out in the regime where the instanton 
density ~ exp{'~^ N) is parametrically suppressed at = oo. From eq. ([6]) 
one can obtain the following expression for instanton density in vicinity of 



where a is 11/3 at one loop level, but the perturbative corrections could be 
large and they may considerably change this numerical coefficient. Such a 
behavior ([9]) does imply that the dilute gas approximation is justified even 
in close vicinity of /Zc as long as H—H^ ^ A., in this case the diluteness 
parameter remains small. We can not rule out, of course, the possibility 
that the perturbative corrections may change our numerical estimate for 
/ic. However, we expect that a qualitative picture of the phase transition 
advocated in this paper remains unaffected as a result of these perturbative 
corrections in dilute gas regime. 

c. Once He is fixed one can compute the entire segment of the phase tran- 
sition line fJ-c{T) for relatively small T. Indeed, in the dilute gas regime 





(9) 
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at /i > /ic the T dependence of the instanton density is determined by a 
simple insertion ~ exp[— 2/3A'^7r^T^p^] in the expression for the density, see 
e.g. the review paper^ . The result in the leading loop order HciT) can be 
presented as follows, 



4. QCD defects via holography 

Our discussions in the previous sections were based on the instanton cal- 
culus in vicinity of the phase transition. The instanton in the holographic 
description can be identified with the DO brane extended along X4,, see 
j.g£gi,3,4 £qj. notations and introduction into the subject. The holographic 
picture gives a very natural way to represent the deconfinement phase tran- 
sition as the Hawking-Page phase transition in which case the two metrics 
with the same asymptotics get interchanged. As our main goal in this talk 
is to understand the phase transition, it would be interesting to study other 
gauge configurations from the holographic viewpoint. In what follows we 
limit ourselves only with two examples. For a more complete list of different 
defects which can be studied by this technique, see the original paper. ^ 
particle 

There are several possible embeddings of D4 branes. One possibility is to 
wrap D4 brane around 5^ and extended along the Euclidean time r, see^ 
for the details. The key point here is that due to this wrapping around 
the "electric charge" N is induced on the D4 brane. The resulting 3 dimen- 
sional system is a static and topologically stable configuration. If we had 
quarks we could make a baryon by attaching quarks to open strings. In the 
pure YM case there are no flavor branes describing quarks to make a gauge 
invariant object. The most simple way to achieve this goal is to add N DO 
branes yielding the D0-D4 open strings. Hence we get the D4 particle which 
is not sensitive to the ^-term. It is important that this configuration is well- 
defined below the critical temperature T < Tc in contrast with instantons 
discussed in^''* when DO branes were sensitive to the 9 term and were well 
defined above the critical temperature T > T^. We can combine the " D4 
particle" with "D4 anti-particle" to form a gauge invariant object, see Fig 
1. The mass of this object scales as ~ TV and is much heavier than the usual 
glueballs. In fact, one can construct gauge invariant objects with any even 
number of vertexes such that the total charge vanishes. It is a new family 
of glueballs with mass ^ NV , where V is the number of D4 and anti-D4 
particles which form a desired configuration. It is amusing that such kind 



Mc(T) = Me(T = 0) 1 



s/NT « iJL^. (10) 



3iV//x2(r = 0)J' 
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^ QCD Junction Buckyballs : (JJ )y/2 

Fig. 1. Buckyballs with N = 3 and V = 8,24,48, 120 from.^^ The construction com- 
bines " D4 particle" with "D4 anti-particle" to form a gauge invariant object with zero 
baryon charge. The mass of this object scales as ~ N. 



of structure in QCD had been previously discussed^^ motivated by the dis- 
covery of the carbonic FuUerenes Ceo and C70 in 1985, which are nano-scale 
objects. ^'^ The QCD objects, similar to the carbonic FuUerenes with femto- 
meter scale were named Buckyballs. It has been also demonstrated that the 
"magic" numbers for Buckyballs are V ~ 8, 24, 48, 120 which correspond to 
the most symmetric, and likely, most stable configurations.^^ The proper- 
ties of this configuration are not sensitive to 6. The possibility to discover 
such kind of configurations at RHIC was discussed in.^^ 

D0-D2 . 

In this subsection we want to address the following question: what happens 
to instantons (represented by DO branes in holographic description) in the 
confined phase? Naively, one could think that the system becomes unstable 
at T < Tc, and therefore the instantons simply disappear from the system. 
However, one should speak about effectively zero action for formation of 
such kind of objects rather than about their instability. Therefore, numerous 
number of these objects can emerge in the system without any suppression. 
In fact, it has been argued in^ that this is precisely what is happening when 
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one crosses the phase transition hne from above. 

It turns out that there are configurations in holographic picture which 
exactly describe dissociation of the instanton into N different objects with 
fractional topological charge 1/A''. The construction is described in^ and 
requires N additional D2 domain walls localized at some points along .T4 co- 
ordinate. In this construction an object with a fractional topological charge 
1/N may emerge. Indeed, one can follow the construction of ref.^^ for SUSY 
case when D2 branes located symmetrically split the instanton into TV 
constituents stretched between pairs of domain walls. Each constituent has 
fractional instanton number 1 /N as well as the fractional monopole number 
and has no reason to condense. In our system we have precisely appropri- 
ate D2 branes which are needed for this construction. These monopoles are 
instantons in the 3d gauge theory on the D2 worldvolume theory which 
involves the scalar corresponding to the position of D2 branes. As we spec- 
ulate below these magnetic monopoles may play an important role in the 
region close to the phase transition < |T — Td/Tg < 1/N. 

5. Conclusion. Speculations. 

A general comment on this proposal can be formulated as follows. Our 

conjecture which relates two apparently unrelated phenomena (phase tran- 
sition vs sharp changes in 9 behavior) implicitly implies that topological 
configurations which are linked to 6 must play a crucial role in the dynam- 
ics of the phase transition. For T > Tc such configurations are well-known: 
they are dilute instantons with density ~ e~''^ cosO. We presented argu- 
ments in^ (see also earlier references therein) suggesting that at T < Tc 
the instantons do not disappear from the system, but rather dissociate into 
fractionally charged constituents, the so-called instanton quarks. In this 
sense the phase transition can be understood as a phase transition between 
molecular phase (deconfined) and plasma phase (confined) of these frac- 
tionally charged constituents. A similar conclusion on sharp changes in 6 
behavior at T = Tc was also observed in ref.-'^^ where the authors studied 
the D2 branes in confined and deconfined phases at T ^ 0. The topological 
objects (sensitive to 6) were identified in ref.^^ as magnetic strings. 

If the picture advocated in the present work about the nature of the 
transition turns out to be correct, it would strongly suggest that fraction- 
ally charged constituents (which carry the magnetic charges as discussed 
in^) may play a very important role in dynamics in deconfined phase in close 
vicinity of the transition < (T— Tc) < For large TV this region shrinks 
to a point, however for finite TV it could be an extended region in tempera- 
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tures. In this region the instantons are not formed yet, and our semiclassical 
analysis is not justified yet as eq. ^ suggests. However, the constituents in 
this region are already not in condensed form. Therefore they may become 
an important magnetic degrees of freedom which may contribute to the 
equation of state similar to analysis on wrapped monopoles in ref.^^ 

To conclude this talk I want to make a few comments on two recent 
papers^^'""^* where the authors advocate the picture similar to the one pre- 
sented here. 

I start with.^^ In that work the authors consider a specifically deformed 
SU (N) gluodynamics at T ^ 0. It has been shown that such a deformation 
supports a reliable analysis in the weak coupling regime in the confining 
phase. The results of the corresponding calculations imply that the rele- 
vant degrees of freedom in the confined phase are the self dual magnetic 
monopoles with action and topological charges Q — ±1/N which are 
precisely the features of the instanton quarks discussed above. 

In contrast, the starting point of ref.^® is semiclassical calculations in the 
background of calorons where the weak coupling regime can not be guaran- 
teed. While the calculations are semiclassical in nature, and therefore, can 
not be trusted in the strong coupling regime, still, the corresponding anal- 
ysis shows how well localized instantons with integer topological charges 
at T > Tc may dissociate into the fractional constituents at T < Tc, and 
become the key players in the confining phase. This is precisely the picture 
we are advocating in the present work based on analysis of sharp 9 changes 
at iV = oo. It is impressive how complicated semiclassical calculations car- 
ried out in^* lead to the expression for the vacuum energy E^ac cos{j^) 
advocated in^^ using completely different technique. 

Our technique does not allow us to make any dynamical calculations 
in this phase as all color degrees of freedom have been integrated out in 
obtaining the low energy effective lagrangian. In other words, we can not 
study the dynamics of fractionally charged constituents in contrast with 
papers^^'^^Q However, the fact that the constituents carry fractional topo- 
logical charge can be recovered in our approach because the color- 
singlet rj' field enters the effective lagrangian as cos(^^) and serves as a 
perfect probe of the topological charges of the constituents. One should also 
emphasize that the procedure of the recovering of the fractional topological 



''In particular, we do not see a beautiful picture of a multi-component color Coulomb 
plasma with nearest-neighbor interactions in the Dynkin space advocated in.^^'^* Still, 
we do see the color- singlet Coulomb interaction of the fractionally charged H/N con- 
stituents due to r]' at very large distances where color already confined. 
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charge 1/N (which has been used here) is not based on the weak coupling 
expansion. 

It was great pleasure to participate in the Workshop honoring 60th 
anniversary of Misha Shifman. Happy birthday, Misha! 

This work was supported, in part, by the Natural Sciences and Engi- 
neering Research Council of Canada. 
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